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ABSTRACT 

GAIA is the “super-Hipparcos” survey satellite selected as a Cornerstone 6 mission 
by the European Space Agency. GAIA can measure microlensing by the brightening 
of source stars. For the broad G band photometer, the all-sky source-averaged photo¬ 
metric optical depth is ~ 10~ 7 . There are ~ 1300 photometric microlensing events for 
which GAIA will measure at least one datapoint on the amplified lightcurve. GAIA 
can also measure microlensing by the small excursions of the light centroid that occur 
during events. The all-sky source-averaged astrometric microlensing optical depth is 
- 2.5 x 10“ 5 . Some r\j 25000 sources will have a significant variation of the centroid 
shift, together with a closest approach, during the lifetime of the mission. This is not 
the actual number of events that can be extracted from the GAIA dataset, as the false 
detection rate has not been assessed. 

A covariance analysis is used to study the propagation of errors and the estimation 
of parameters from realistic sampling of the GAIA datastream of transits in the along- 
scan direction during microlensing events. The mass of the lens can be calculated to 
good accuracy if the lens is nearby so that angular Einstein radius Oe is large; if the 
Einstein radius projected onto the observer plane fg is about an astronomical unit; if 
the duration of the astrometric event is long (> 1 year) or if the source star is bright. 
Monte Carlo simulations are used to study the ~ 2500 events for which the mass can 
be recovered with an error of < 50%. These high quality events are dominated by disk 
lenses within a few tens of parsecs and source stars within a few hundred parsecs. 
We show that the local mass function can be recovered from the high quality sample 
to good accuracy. GAIA is the first instrument with the capabilities of measuring 
the mass locally in very faint objects like black holes and very cool white and brown 
dwarfs. 

For only ~ 5% of all astrometric events will GAIA record even one photometric 
datapoint. There is a need for a dedicated telescope that densely samples the Galactic 
Centre and spiral arms, as this can improve the accuracy of parameter estimation by 
a factor of ~ 10. The total number of sources that have an astrometric microlensing 
variation exceeding the mission target accuracy is ~ 10 5 . The positional measurement 
of one source in every twenty thousand is affected by microlensing noise at any instant. 
We show that microlensing is negligible as an unbiased random error source for GAIA. 

Key words: gravitational lensing - astrometry - Galaxy: stellar content - Galaxy: 
structure - dark matter 


1 INTRODUCTION 


precedented precision on over a billion objects in our Galaxy 
alone. 


GAIA is the European Space Agency (ESA) satellite now 
selected as a Cornerstone 6 mission as part of the Science 
Program It is a survey satellite that provides multi-colour, 
multi-epoch photometry, astrometry and spectroscopy on all 
objects brighter then V ss 20 (e.g., ESA 2000; Perryman 
et al. 2001). The dataset is huge with information of un- 


GAIA is the successor to the pioneering Hipparcos satel¬ 
lite, which flew from 1989 to 1993. The Hipparcos program 
was both simpler and smaller: it measured only 10 5 rather 
than 10 9 objects and it was provided with a target list rather 
than being left to determine its targets for itself, as GAIA 
will. The Hipparcos final results were released only when the 
mission was complete, whereas many of the science goals of 
GAIA, especially for bursting or time-varying phenomena 
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mas 


Figure 1. This shows the relative right ascension and declination 
of a source with (solid line) and without (dashed line) a microlens- 
ing event. The trajectory is shown over the GAIA mission lifetime 
of 5 years. Note that the deviations caused by microlensing are 
present well before and after the time of maximum of the event. 
(The lens is at 150 pc, the source at 1.5 kpc, the transverse ve¬ 
locity is 70kms~ 1 and the impact parameter u is 1.5. The lens 
has mass 0.5Mg.) 


like microlensing, may require an early analysis and release 
of preliminary data. 

GAIA carries out continuous scanning of the sky. The 
satellite rotates slowly on its spin axis, which itself precesses 
at a fixed angle to the Sun of 55°. As GAIA rotates, light 
enters the entrance chambers, reflects off mirrors and falls 
on the focal planes of one of three telescopes, two of which 
(ASTRO-1, ASTRO-2) measure the positions of stars, one of 
which (SPECTRO) performs spectroscopy. GAIA observes 
in three directions along a great circle simultaneously. The 
astrometric parameters of stars are recovered from the time 
series of the one-dimensional transits distributed over the 
five year mission lifetime. 

A small fraction of the objects monitored by GAIA will 
show evidence of microlensing. GAIA can observe microlens¬ 
ing by measuring the photometric amplification of a source 
star when a lens and a source are aligned. This is the ap¬ 
proach followed by the large ground-based microlensing sur¬ 
veys like MACHO, EROS, OGLE and POINT-AGAPE (see 
e.g., Alcock et al. 1997; Aubourg et al. 1995; Udalski et 
al. 1994, Auriere et al. 2001). GAIA is inefficient at discov¬ 
ering photometric microlensing events, as the sampling of 
individual objects is relatively sparse (there are a cluster of 
observations once every two months on average). 

However, there is a much more powerful strategy avail¬ 
able to GAIA. Although the two images of a microlensed 
source are unresolvable, GAIA can measure the small devi¬ 
ation (of the order of a fraction of a milliarcsec) of the cen¬ 
troid of the two images around the trajectory of the source. 
Astrometric microlensing is the name given to this excur¬ 
sion of the image centroid. The cross section of a lens is 
proportional to the area it sweeps out on the sky, and so 
to the product of lens proper motion and angular Einstein 
radius. Each of these varies as the inverse square root of 
lens distance, so the signal is dominated by nearby lenses. 
The detection of astrometric microlensing events for pointed 
observations has been considered many times before (e.g., 
Walker 1995; Miralda-Escude 1996; Paczynski 1996; Boden, 
Shao & van Buren 1998). Once alerted by a ground-based 
survey, the events are followed by narrow angle differential 


astrometry. The detection of astrometric microlensing using 
observations from a scanning satellite has been considered 
before by Hog, Novikov & Polnarev (1995) in the context of 
the proposed ROEMER mission (a forerunner of GAIA). 

We will show that the all-sky source-averaged astromet¬ 
ric microlensing optical depth is ~ 2.5 x 10 —5 , over an or¬ 
der of magnitude greater than the photometric microlensing 
optical depth. There are two main difficulties facing GAIA 
in exploiting this comparatively high probability. First, the 
astrometric accuracy of a single measurement by GAIA de¬ 
pends on the source magnitude and quickly degrades at mag¬ 
nitudes fainter than G ~ 15. Second, GAIA provides a time- 
series of one-dimensional astrometry. The observed quantity 
is the CCD transit time for the coordinate along the scan. 
This is the same way the Hipparcos satellite worked (see 
The Hipparcos and Tycho Catalogues, ESA (1997), volume 
3, section 16). From the sequence of these one-dimensional 
measurements, the astrometric path of the source, together 
with any additional deflection caused by microlensing, must 
be recovered. 

This paper assesses the microlensing signal that will be 
seen by GAIA. Section 2 gives the formulae for the photo¬ 
metric and astrometric microlensing optical depths. These 
are used to estimate the total number of microlensing events 
that GAIA can measure. In Section 3, a covariance analysis 
is used to demonstrate the effects of error propagation on 
the recovery of the parameters of microlensing events. For 
the purposes of GAIA, we show that the disk stars within 
a few hundred parsecs of the Sun are the most important 
source and lens populations. Section 4 describes the extinc¬ 
tion law and the Galaxy model used to generate microlens¬ 
ing events for our Monte Carlo simulations. The synthetic 
data are sampled with GAIA’s scanning law and realistic 
errors are applied to provide the one-dimensional astromet¬ 
ric datastreams. Section 5 discusses the results of the sim¬ 
ulations, both for the entire sample of events and for the 
subset of gold-plated events whose parameters can be re¬ 
covered to good accuracy. Finally, Section 6 examines the 
overall strategy for identification of events, which need to 
be distinguished from other forms of astrometric deviation 
such as binary companions. 

2 ASTROMETRIC MICROLENSING 
2.1 General Principles 

Let us consider the lensing of a luminous source by a dark 
point lens. (Of course, the lens may be luminous, but we de¬ 
fer consideration of this till Section 6). The angular position 
of the source 0 B can be written as 

@s(t ) — @5,0 + p s £ + Ps, ( 1 ) 

where @ B ,o is the zero-point, p a is the source proper motion 
and P s is the source parallax. Similarly, the angular separa¬ 
tion between source and lens is given by 

@si(t) = @s\,o + M s it + Psi, (2) 

where the proper motion p sl and parallax P a i are of the 
source relative to the lens. 

Microlensing induces an additional shift of the source 
centroid. Letting S(f) be the image centroid, we have (e.g., 
Walker 1995; Dominik & Sahu 2000) 
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G (ill mag) 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

<7 p (in mmag) 

9 

9 

9 

9 

9 

10 

12 

16 

23 

36 

60 

cr a (in pas) 

30 

30 

30 

40 

60 

90 

150 

230 

390 

700 

1400 


Table 1. This table lists the mean accuracy in photometry cr p and in position <r a versus G band magnitude for the GAIA satellite. Note 
that cj a is the accuracy of a single astrometric measurement, not the target accuracy at the end of the GAIA mission, which is better by 
a factor of ~ 10. The values are approximate sky averages, adapted from Tables 7.3 and 8.2 of ESA (2000). 



Figure 2. Upper panel: Astrometric shift of the microlensing 
event of Figure 1, as seen by a barycentric (dotted line) and a 
terrestrial observer (solid line). Lower panel: Simulated data in¬ 
corporating typical sampling and astrometric errors for GAIA. 
Also shown for comparison are the theoretical trajectories of the 
source with (grey line) and without (dashed line) the event. The 
insets show the deviations at the beginning and the midpoint of 
this high signal-to-noise event. (The accuracy cr a of the astrome¬ 
try is 300 //as, corresponding roughly to a 17th magnitude star). 


= ° = w "w + r (3) 

Here, the angular Einstein radius 6e is related to the lens 
mass M by 


Figure 3. Same as Figure 2, but the lens distance is 1 kpc. This 
shows a microlensing event with a 5v^2<Va variation of the centroid 
shift. GAIA will be able to measure the astrometric deviation 
but will not be able to extract any useful information on the 
microlensing parameters. 


The characteristic lengthscale in microlensing is the 
Einstein radius Re = A$e- It is the linear size of the angu¬ 
lar Einstein radius in the lens plane. For applications in the 
Galaxy, we find that 


Re = 9 au 



A 

10 kpc 



A 

DA 


(5) 



where 7r s i = |P a i|- Figure |l] shows the astrometric trajectory 
of the same source with and without a microlensing event. 
In this case, the lens is just 150 pc away from the observer, 
while the source is a disk star 1.5 kpc away. The microlensing 
event causes the additional excursion superimposed on the 
parallactic and proper motion of the source. 


and so Re is of the typical size of a few astronomical units 
(aus). Here, D\ and D s are the distances from the observer to 
the lens and source respectively. The characteristic timescale 
in photometric microlensing is the Einstein crossing time Ie, 
which is the time taken for the source to cross the Einstein 
radius. 
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Figure 4. Upper panel: An all-sky map of the source-averaged as¬ 
trometric optical depth. Meridians of galactic latitude are shown 
at 60° intervals, parallels of longitude at 30° intervals. The con¬ 
tour levels are marked. The number within each contour refers 
to the mean value of the optical depth inside the shaded region 
enclosed between the contours. All numbers are in units of 10~~ 5 . 
Middle panel: As the upper panel, but for the photometric mi- 
crolensing optical depth. All numbers are in units of 10“ 7 . Lower 
panel: An all-sky map of the star counts. The number within each 
contour refers to the total number of stars inside the shaded re¬ 
gion enclosed between the contours. To find the instantaneous 
number of events in any region, we multiply the number of stars 
by the average optical depth. (For details of the source and lens 
populations and extinction model, see Section 4). 


2.2 The Centroid Shift 


The centroid shift points away from the dark lens as seen by 
the source. The dimensionless function 


u = 


flsl 


( 6 ) 


is the angular separation in units of the angular Einstein 
radius. For u —> oo, the centroid shift falls off like 


S=|S|~-. (7) 

u 

This can be compared with the photometric magnification 
A, which falls off like 

A ~ 1 + ^ < 8 > 

These asymptotic results illustrate one of the important dif¬ 
ferences between astrometric and photometric microlensing. 
The centroid shift falls off much more slowly than the mag¬ 
nification, so that the cross-section for astrometric events 
is much larger than for photometric events (e.g., Paczynski 
1996; Miralda-Escude 1996). 

In the absence of the proper and parallactic motion, the 
absolute value of the centroid shift is 


e = 


V u o + i 2 p 

—---C/e 

Uq + t 2 + 2 


(9) 


Here, uo is the value of the dimensionless distance at the 
time of closest approach to and 


_ t — to _ — to) 

tE 0 e 


( 10 ) 


The shift can be decomposed into components parallel and 
perpendicular to the direction of motion of the lens relative 
to the source, namely 


9 II = 


Ug+t 2 + 2 


9e, 


e± = 


uo 


9e- 


Up + 1 2 T 2 


( 11 ) 


As pointed out by Walker (1995), the centroid shift as seen 
by a barycentric observer is an ellipse. For u —> 0, the ellipse 
becomes a straight line, while for u —> oo, it becomes a circle. 

Figure)^ shows two further views of the same microlens¬ 
ing event presented in Figure [lj The upper panel shows the 
right ascension and declination recorded by a barycentric 
and a terrestrial observer (or equivalently a satellite at the 
L 2 Lagrange point, like GAIA). The proper and parallactic 
motion of the source have been subtracted out. However, 
GAIA does not provide the data in such a clean form, as 
it really measures a series of one-dimensional transits of the 
two-dimensional astrometric curve. The lower panel shows 
the event as seen by GAIA. The simulated datapoints have 
been produced by generating random transit angles, and 
sampling the astrometric curve according to GAIA’s scan¬ 
ning law for the ASTRO-1 and ASTRO-2 telescopes (using 
programs freely available from L. Lindegren at Lund Obser¬ 
vatory). The transits are strongly clustered, as GAIA spins 
on its axis once every 3 hours and so may scan the same 
patch of sky four or five times a day. The transit angle is 
the same for all transits in such a cluster, but changes ran¬ 
domly from cluster to cluster. Gaussian astrometry errors 
with standard deviation <r a = 300 ^tas have been added to 
the simulated datapoints. The two insets show the astromet¬ 
ric deviations at the beginning and at the maximum of the 
event, from which it is clear that GAIA can detect that a 
microlensing event has occurred. Figure shows the same 
microlensing event, except that the lens distance has been 
increased. This causes the relative parallax 7r a i, and conse¬ 
quently the angular Einstein radius #e, to diminish. The 
upper panel of Figure |j| shows that the parallactic deviation 
of the curve is also smaller than before. It is already clear 
from the lower panel that it will be difficult to extract phys- 


© 0000 RAS, MNRAS 000, 000-000 



















































































Astrometric Microlensing with the GAIA satellite 5 


Number of transits in Astro—1,2 






Ecliptic latitude, b 


Figure 5. This shows the number of transits made by the 
ASTRO-1 and ASTRO-2 telescopes as a function of ecliptic lat¬ 
itude. The figure is drawn by choosing 5000 random directions 
on the sky, and computing the ecliptic latitude and number of 
transits using software supplied by L. Lindegren. It should be re¬ 
membered that the number of transits may be misleading, as the 
transits are strongly clustered into groups of between two and 
five. 


ical information from some of the events, particularly those 
with smaller Einstein radius. 


2.3 Optical Depths 


The probability r of detecting a microlensing event with 
certain characteristics is just 

r = D s J^ dxJ (12) 

where E is the area in the lens plane for which the projected 
source positions yield the required effect and f(M) is the 
mass function or the number of lenses with mass in the range 
M to M + dM. 

GAIA measures displacements only along the scan, so 
we incorporate a factor of y/2 to take us from GAIA’s one¬ 
dimensional astrometric accuracies to the two-dimensional 
accuracies. The centroid shift varies by more than 5v2cr a 
when the projected position of the source lies within a circle 
of radius u^Re, where (Dominik & Sahu 2000) 


TjifeV 

5V2a a D\ 


Tiife fe 

5\/2<7- a tE 


(13) 


Here, Tiif e is the lifetime of the satellite and v is the trans¬ 
verse source velocity in the lens plane. Note that 5y/2a a is 
enough for GAIA to make a convincing measurement of the 
shift, but not usually enough for detection of a microlensing 
event. 

Instead, the minimum criterion needed to detect the 
event in the GAIA database is a variation of centroid shift 
larger than 5\/2<r a together with a closest approach between 
lens and source happening with a time Tiif e . The area in the 
lens plane giving rise to events with these characteristics is 
E = 2u a R E Ti i{e v. Using (^), the probability is (Dominik & 
Sahu 2000) 



where 


(M~ 1/2 ) = / M~ 1/2 f(M)dM. (15) 

Jo 

We refer to the quantity ([ifj) as the astrometric microlensing 
optical depth. 

By contrast, the probability of a microlensing event 
with an amplification greater than that corresponding to 
the threshold of photometric accuracy 5a p is (e.g., Griest 
1991, Dominik & Sahu 2000), namely 


4nG 


u 2 D 2 / dx p(x)x{ 1 — x), 


where 


2 

u„ = 


P y/TT4 


-2 

P 


- 2 . 


(16) 


(17) 


Here, A p is the magnification corresponding to the threshold 
photometric accuracy 5<r p . We refer to the quantity (|l^) as 
the photometric microlensing optical depth. For GAIA, both 
the astrometric and photometric accuracies depend on the 
magnitude of the source as listed in Table 

The Einstein crossing time t E is the time taken for the 
source to cross the Einstein radius. This is not the dura¬ 
tion of a photometric event, which is the time taken for the 
source to cross a circle of radius u p Re in the lens plane. 
The duration of an astrometric event is the time for which 
the centroid shift is greater than the threshold 5\Z2a a . Do¬ 
minik & Sahu (2000) show that this is approximately the 
time taken to cross a circle of radius I?E#E/(5\/2cr a ). So, the 
duration of an astrometric event is 


t E 


9e 

5\/2cr a ' 


(18) 


For GAIA, this is typically a factor of two times longer than 
the duration of a photometric event. 

Figure ^ show contours of astrometric and photomet¬ 
ric microlensing optical depth together with starcounts. The 
maps assume a standard model for the sources and lenses in 
the Galaxy, together with a luminosity function and extinc¬ 
tion law (described in Section 4). They have been produced 
assuming that the relative source-lens velocity in the lens 
plane is ~ 140kms -1 (see Figure |l^). Extinction is an im¬ 
portant effect for GAIA’s microlensing capabilities, as the 
accuracy of both the astrometry and photometry depends 
on source magnitude. Taking into account the effects of ex¬ 
tinction, the all-sky averaged value of the astrometric opti¬ 
cal depth is 2.5 x 10~ 5 . Here, the averaging is performed by 
weighting the optical depth with the starcount density. Re¬ 
gions like the Galactic bulge are very heavily weighted, and 
so the average astrometric optical depth is of the same order 
of magnitude as the typical optical depth towards the cen¬ 
tral part of the Galaxy. There are ~ 10 9 stars brighter than 
V = 20 in our model; the same is true for the Galaxy (Mi- 
halas & Binney 1981). This means that, during the GAIA 
mission, there are ~ 25000 astrometric microlensing events, 
which have a variation of the centroid shift greater than 
5y/2 cr a together with a closest approach during the mission 
lifetime. This number can be read off Figure by multiply¬ 
ing the mean astrometric optical depth within the central 
contour of the upper panel by the total number of stars 
within the central contour of the lower panel. Some of these 
displacements cannot be identified by GAIA as microlens¬ 
ing events, first because the signal-to-noise will often be low 
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and second because any identification algorithm must not 
generate too many false detections. The middle panel of 
Figure ^Jshow contours of photometric microlensing optical 
depth again including the effects of extinction. The all-sky 
averaged photometric optical depth is ~ 1.2 x 1CP 7 . Let us 
recall that, from eqs ( |Ig| - |I7| ), a photometric event occurs 
whenever there is 5 ct p variation. The typical event duration 
is ~ 2 months, so that there are a total of ~ 3600 photomet¬ 
ric microlensing events during the GAIA mission. However, a 
very substantial number of these events will be undetectable. 
GAIA’s sampling is sparse compared to ground-based pro¬ 
grams like MACHO, EROS or POINT-AGAPE, so many of 
these events will be missed. 


3 ERROR PROPAGATION 
3.1 Covariance Analysis 

Our first task is to understand the propagation of errors, 
so that we can assess how many of the ~ 25000 potential 
astrometric microlensing events are useful. 

Since most events will not be detected photometrically 
by GAIA, it will not usually be possible to combine pho¬ 
tometric together with astrometric data to analyze a mi¬ 
crolensing event. The measured quantity that will generally 
be provided by GAIA is the source displacement along the 
scan. This contains information on the microlensing event, 
but is contaminated with the source parallactic and proper 
motion as well. (We assume that the GAIA datastream has 
already been corrected for aberration due to satellite mo¬ 
tion and gravitational deflection caused by Solar System 
objects). 

The components of S resolved with respect to right as¬ 
cension a and declination <5 are 

(s a ,s t ) = e s + ^— 0E. (19) 

Ip + 2 

The angular position of the source 9 S contains the contribu¬ 
tions from the source proper motion fi s and source parallax 
P s . We note that 

( t cos <j> — uo sin <j> — n s iP a 8 ^, 1 
t sin <f> + uo cos </> — 7T s i PsS^ 1 

where <j> is the proper motion angle (Gould & Salim 1999). 
We have also expanded the parallaxes P s and P s i as 

Ps = 7T S (Pot, Ps) , Psl = TTsl [Pa, Pi) , (21) 

so that 7r s and 7r s i are the absolute values and P a and Ps are 
the direction cosines (given in van der Kamp 1967). 

The directions along and perpendicular to the scan are 
related to right ascension and declination by a rotation, 
which depends on where the satellite is pointing. As GAIA 
continuously scans the sky, this angle can take any value 
for the same source. So, formulae (pli|)-(|2lj) hold good in the 
scan coordinate system as well, after rotation by a random 
angle to take us from the right ascension and declination to 
the along-scan and across-scan coordinates. GAIA measures 
positions in the along-scan direction only. 

We use this microlensing model to create synthetic as¬ 
trometric measurement sets. There are a total of 11 param¬ 
eters which we wish to compute from the data, namely: /x s , 



Parallax uncertainty vs. projected Einstein radius 



Figure 6. This shows the percentage error in estimation of the 
relative parallax as a function of the projected Einstein radius 
rg for different astrometric event durations and a given accuracy 
cr a = 150 /ms. The events have varying source and lens distance, 
as well as the transverse velocity. (All the other parameters are 
as in Figure 1, except u = 1). 


Parallax uncertainty vs. accuracy 



^as 


Figure 7. This shows the percentage error in estimation of the 
relative parallax as a function of astrometric accuracy for different 
angular Einstein radii (?e • We vary the Einstein radius by varying 
the source and lens distance at fixed transverse velocity. (All the 
other parameters are as in Figure 1, except u = 5). 


7t s , Uo, to, t e, 4>, 9e, 7t s i and 0 s ,o- To study parameter estima¬ 
tion performance, we evaluate the covariance matrix dj (see 
e.g., Boutreux & Gould 1996; Gould & Salim 1999). This is 
defined as the matrix C = B _1 , where 



d0 a (tk) d0 a (tk) 
dai daj 


( 22 ) 


Diagonal elements of the matrix ^fcii are the individual un¬ 
certainties in the parameters assuming the standard 
fitting procedure. Here, fli,...,on are the 11 parameters 
stated above and er a is the astrometric accuracy for a given 
star, which depends on the source magnitude and is recorded 
in Table [jj 

Note that tk are the times at which GAIA samples the 
one-dimensional astrometry datastream. In our calculations, 
these times are always generated from GAIA’s scanning law 
using Lindegren’s software. GAIA’s scanning law is quite 
complicated, as shown in Figure |E|. The number of times 
GAIA samples an object N is between 100 and 450, depend¬ 
ing on its ecliptic latitude (Lindegren 1998). The number of 
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Uncertainties vs. impact parameter 



Figure 8. This shows the percentage error in estimation of the 
microlensing parameters as a function of impact parameter u. Ac¬ 
curate recovery of the relative parallax 7r s j and the proper motion 
angle <j> is harder, whereas recovery of the angular Einstein radius 
Oy and the Einstein crossing time t E is easier. (The event has the 
same parameters as Figure 1 with <r a = 150 //as.) 


Uncertainties vs. sampling 



Figure 9. This shows the percentage error as a function of num¬ 
ber of datapoints for a typical microlensing event. Provided there 
are > 150 measurements, then the sampling is not the critical 
factor. (The event has the same parameters as Figure 1, but 
(T a = 100 lias). 

transits is greatest, and so the final astrometry is most ac¬ 
curate, for objects close to an ecliptic latitude ~ 35°. For 
a typical location, the number of transits is between 100 
and 200. These transits are usually clustered into groups of 
between 2 and 5, so the average effective number of inde¬ 
pendent datapoints is ~ 40. 

3.2 Errors in Estimates of the Astrometric 
Microlensing Parameters 

Although the relative parallax 7r s i is quite hard to estimate, 
it is valuable since it provides a unique solution for the mass 
of the lens M when combined with # E (see eq. (Q)). If in addi¬ 
tion the source distance D s can be estimated for example, 
if GAIA itself measures the source parallax - then we have a 
complete solution for all the microlensing parameters. Let us 
recall that in the ground-based MACHO, EROS and OGLE 
programs, there is not a single microlensing event for which 


the parameters can be derived unambiguously from the data 
alone. Consequently, the location of the lenses responsible 
for the events towards the Magellanic Clouds is a matter of 
considerable controversy (e.g., Sahu 1994, Gould 1995, Al- 
cock et al. 1997, Evans & Kerins 2000, Zhao & Evans 2000). 
GAIA offers the promise of providing a sample of lensing 
events for which all the parameters can be extracted with¬ 
out any modelling assumptions. However, some caution is 
needed, as we already known in the context of the Space In¬ 
terferometry Mission (SIM) by Gould & Salim (1999) that 
astrometric measurements alone do not necessarily guaran¬ 
tee accurate estimation of the relative parallax. 

For which events can we recover the relative parallax to 
good accuracy from the GAIA data? First of all, an event 
must last an adequate time for the microlensing shift to be 
distorted by parallactic movement of the lens. Thus, the er¬ 
ror will depend on the duration of the astrometric event. Sec¬ 
ond, the amplitude of the distortion is dictated by the Ein¬ 
stein radius projected onto the observer plane (e.g., Gould 
& Salim 1999, Gould 2000), namely 

r E = — • (23) 

TTsl 

For close lenses (Di ss 100 pc), it follows that 1—D\/D S —♦ 1, 
and so the projected Einstein radius becomes 

™ (24) 

For a measurable distortion, we require r E to be about an as¬ 
tronomical unit or smaller. If it is too large, then the Earth’s 
motion about the Sun has a negligible effect. Accuracy in the 
relative parallax is a trade-off between duration and lens dis¬ 
tance. 

Figure ^ shows the uncertainty in 7r s i as a function of 
r E for different event durations t ae , while the companion 
Figure Q shows the uncertainty in 7r s i as a function of astro- 
metric accuracy for different # E . It aids accurate recovery of 
the relative parallax if the angular Einstein radius is large. 
The error in 7r s i degrades for distant lenses. The degrada¬ 
tion is worse for the shorter duration events. Note that the 
close lenses {D\ « lOOpc) typically have a crossing time 
t E « lau/150 knis' 1 « 10 days, and so the astrometric 
event duration is typically 200 days. It is the close lenses 
with longer timescales that provide the most propitious cir¬ 
cumstances for measuring 7r s i from the data. 

Figure |i] shows the uncertainties in # E , t E , 7r s i and (f> as 
a function of the dimensionless impact parameter. A typi¬ 
cal value of u for GAIA is 11 or so. For large u, the par¬ 
allactic distortion is superposed upon a small eccentricity 
ellipse (see eq. ©)• For small u, the eccentricity of the 
ellipse becomes large and the parallactic distortion makes 
the curve thinner and more elongated. This can be an al¬ 
most degenerate situation for estimating the microlensing 
parameters from one-dimensional astrometry. The percent¬ 
age error in (j> and 7r s i is larger than the error in # E and 
t E . This is because 7r s i controls the shape of the astrometric 
curve, <j> controls the orientation, while # E and t E the size. 
Given one-dimensional astrometry, it is more difficult to re¬ 
cover the details of the shape than the size. The rise and 
fall in error at u « 3 in Figure ^ is a consequence of the 
degeneracy of the one-dimensional astrometry. Asymptoti¬ 
cally, the percentage error always increases with increasing 
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impact parameter. However, there is a regime at small im¬ 
pact parameter when the converse is true. The reason for 
this is as follows. For very close lenses, the parallactic dis¬ 
tortion can become so large that the direction of traversal 
of the astrometric curve is reversed. This degeneracy is in¬ 
herent in microlensing and it is not a consequence of the 
limited information that can be extracted from measure¬ 
ments. The astrometric curve can now correspond to events 
with positive impact parameter and large parallactic effect, 
or negative impact parameter and small parallactic effect. 
There is no way to tell the difference. 

Finally, Figure shows the effects of sampling. The er¬ 
ror in selected microlensing parameters is plotted against 
the number of datapoints for a typical microlensing event. 
In this Figure only, the sampling is not performed according 
to GAIA’s sampling law, but the datapoints are chosen uni¬ 
formly over the mission lifetime, as we wish to illustrate the 
importance of sampling. In fact, GAIA samples every object 
between 100 and 450 times depending on ecliptic latitude. 
So, at first glance, Figure seems to offer reassurance that 
sampling is not a critical factor, as it shows that the error in¬ 
creases swiftly only once the effective number of datapoints 
falls below 40. However, this is misleading because GAIA’s 
datapoints are clustered into groups of between 2 and 5, or 
typically 4. This is disadvantageous, as it reduces the effec¬ 
tive number of datapoints by ~ 4, although the accuracy of 
the astrometry is improved by ~ 2. Figure tells us that 
if the number of independent datapoints falls below ~ 40, 
then sampling is the limiting factor in the extraction of use¬ 
ful parameters. (Our calculations show that the form of the 
figure is only weakly dependent on er a ). In other words, as 
~ 40 is roughly the average number of effective datapoints, 
the sampling is responsible for roughly tripling or quadru¬ 
pling the errors in the estimated parameters for about half 
the dataset. 


Luminosity Functions 



M v 

Figure 10. This shows three luminosity functions, namely the 
Bahcall-Soneira (dashed line), the Reid-Hawley (dotted line) and 
the flat model (full line). The flat model is the standard one in 
our simulations. 


Mass Functions 



l°g(M/M 0 ) 

Figure 11. This shows three mass functions that are flat, rising 
and falling below 0.5 Mq. Above 0.5 Mq, the mass function is de¬ 
rived from the Reid-Hawley LF. Below 0.5 Mq, the mass function 
cannot be derived with certainty from the observations. 


3.3 Errors in Estimates of Masses of Lenses 


We are primarily interested in finding the masses of the 
lenses. The error in the mass is related to the errors in the 
angular Einstein radius and the relative parallax by 



i_ TTsI A /'~*//3 _ \ ^sl 

H-4C(6 l B,7r s i)—-, 

7Tsl f'E TTsl 


(25) 


where C is the correlation coefficient between 0 b and 7r s i. 
Let us recall that the elements of the correlation matrix are 
customarily defined as 





(26) 


The cross-term is important because errors in 0 b and 7r s i are 
strongly correlated for most of the events. This is a small 
term for low impact parameter events, but it becomes im¬ 
portant once u a > 2. The cross-term always decreases the 
uncertainty in the mass. This is why Boden et al. (1998) 
found that eq. (14) of their paper (which neglects the cor¬ 
relation) tends to overestimate the mass error. We always 
use the full formula ( pB| ) to compute the uncertainty in the 
mass in our simulations. 


4 A MODEL OF THE GALAXY 

The goal of the analysis of the microlensing events seen by 
GAIA is to infer properties of the lenses. The likely results 
can best be judged by Monte Carlo simulations. Before this, 
we need to develop a model of the lens and source popula¬ 
tions in the Galaxy. This section describes the ingredients 
of the model in turn. 


4.1 Luminosity Function of Sources 

Let the mass density of sources be p and the luminosity 
function (LF) be <j>. In an element of solid angle dfl, the 
number of stars in a direction (£, b) with distances between 
D s and D s + dD s and with magnitudes between m g and 
m-G + drriG is 

dN = —<j>(M G )dm G p(£, b, D s )D^dD s dQ. (27) 

Po 

Here, £ and b are Galactic coordinates, so that the solid 
angle dfl is just cos b dldb and po is the local mass density. 
As usual, absolute magnitude Mg is related to apparent 
magnitude mo and extinction via 

Mg = m G — 51ogD s — A(D S ,£, b) + 5, (28) 
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Number of stars 
(extinction) 

Number of stars 
(no extinction) 

Luminosity Function 

0.9 x 10 9 

5.0 x 10 9 

flat model 

0.9 x 10 9 

5.0 x 10 9 

Reid, Hawley (2000) 

1.1 x 10 9 

6.1 x 10 9 

Bahcall, Soneira (1980) 


Table 2. The total number of stars with 10 < G < 20 for each 
luminosity function, including and excluding the effects of extinc¬ 
tion. 

where A(D s ,l,b) is the extinction law. Equation ( p?| ) can be 
interpreted as a probability density for the source parame¬ 
ters D s ,£,b and Mg■ 

Three source LFs are shown in Figure [iT| . These are 
valid for the V band, whereas GAIA’s photometry is most 
accurate in the broad G band. As our aim is to provide ap¬ 
proximate numbers in our calculations, the conversion from 
V to G band in the LF will not have a noticeable effect on 
the results. 

The Bahcall-Soneira (1980) LF is a reasonably accu¬ 
rate predictor of the numbers of stars in the Galaxy fainter 
than G ~ 15. The Reid-Hawley LF (2000) is derived from 
a combination of sources, including the Third Catalogue of 
Nearby Stars, the 8 Parsec Sample, as well as Hipparcos 
data. It is the local LF for the stellar disk, but it may suf¬ 
fer from incompleteness at magnitudes fainter than G ~ 16. 
Accordingly, our standard assumption is that the LF is flat 
in the magnitude range 10 < Mg < 20 and follows the Reid- 
Hawley LF when Mg < 10. In this case, the normalisation 
is chosen to be the maximum value of the stellar luminos¬ 
ity function of the Galactic disk as recorded in Table 7.3 of 
Reid & Hawley (2000). The total numbers of stars in the 
range 10 < G < 20 in the three models are given in Table 
The microlensing results reported in this paper are always 
calculated assuming our standard LF. This table tells us the 
rough scaling corrections that we need to convert the results 
to other LFs. 


4.2 Extinction 


the molecular ring at ~ 4 kpc. Nonetheless, it will be good 
enough for our study, as the sources are generally nearby. 

In the Monte Carlo simulations, once the source dis¬ 
tance, direction and absolute magnitude are chosen, then 
the extinction law is used to work out the apparent magni¬ 
tude. We calculate the astrometric accuracy cr a for a given 
source magnitude using the accuracy-magnitude relation for 
GAIA, listed in Table |l]. 

4.3 Density of Sources 

The density of sources at a distance D s in the direction of 
Galactic longitude l and latitude b is 

p = pd + pb- (30) 

Here, pa is the density distribution of the Galactic disk and 
Pb is the density of the Galactic bulge. We use an exponential 
disk model with scalelength Rd ~ 3 kpc and scaleheight 
Zh ~ 350 pc (Gould, Bahcall & Flynn 1997), namely 

pa = poexp (-(R - R 0 )/Rd - \z\/zh) , (31) 

where R is Galactocentric distance, z is the height above or 
below the Galactic plane, Ro is the radial position of the Sun 
and po is the local density of disk stars, which is taken as 
O.IOMqpc -3 (Holmberg & Flynn 2000). It is straightforward 
to establish that 

R 2 = R^+D 2 cos 2 b— 2D s Ro cosbcos£, z = L> s sinb. (32) 

We use the oblate axisymmetric bulge model provided by 
Kent (1992) 

Pb(R,z) = 3.53 K 0 ^ 66 y pc ^ M 0 pc~ 3 , (33) 

where s 4 = R 4 + ( 2 /O. 6 I ) 4 and Ko is a modified Bessel 
function. Kent’s model is a reasonable fit to the Spacelab 
infrared data. However, it is known to be wrong in detail 
as the Galactic bulge is really a triaxial bar (e.g., Binney et 
al. 1991; Hafner et al. 2000). This means that our model - 
like any axisymmetric model - will underestimate the pho¬ 
tometric optical depth as compared to the observations in 
windows close to the Galactic Centre typically by a factor 2 
or so (see e.g., Evans 1994, 1995; Binney 2000). 


We use the standard extinction law 

ro, 


I >50° 


A(D s ,£,b) = 


0.165(1.192-1 tan b\) 

| sin b\ X 

[l-exp(-^s^M)] , 10° < |b| < 50 c 


(29) 


, 7-Ds, 


0 ° < 161 < 10 °. 


For | 6 | > 10°, this is the Sandage absorption law (Sandage 
1972, Chen et al. 1998). The constant ho is the characteristic 
height of the extinction structures, which is ~ 120 pc. For 
\b\ < 10°, a constant differential extinction 7 ~ 0.5 (in mag¬ 
nitudes per kpc) is assumed. This is perhaps a little low, but 
we have chosen it to ensure that the total number of stars 
with magnitudes satisfying 10 < G < 20 is ~ 10 9 , for which 
there is good evidence from starcounts (Table 4.2 of Mihalas 
& Binney 1981). Our extinction model is a reasonably accu¬ 
rate local description, but it does not include features like 


4.4 Mass Function and Density of Lenses 

It is nearby objects that make the largest contribution to the 
astrometric lensing signal. Accordingly, in our Monte Carlo 
simulations, only the disk stars act as lenses. Halo popula¬ 
tions of dark objects like Machos are not included. So, the 
density of lenses is the double exponential disk prescribed 
by eq (^). We generate the masses of the lenses using the 
composite mass function shown in Figure [n|. Above 0.5 Mq, 
we use a mass function (MF) derived from the Reid-Hawley 
LF via 

/(M)=0(M G (M))^, (34) 

Here, Mg(M) is the relationship between mass and absolute 
magnitude, as recorded in Kroupa, Tout & Gilmore (1990). 
This ensures that the high-mass end of the MF (which gives 
most of the light) is compatible with the local LF. The low- 
mass end of the MF is not at all certain and its derivation 
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Parameter 

Whole Sample 

High Quality Sample 

t ae (days) 

32 (163) 

342 (814) 

$E (mas) 

7.9 (20.0) 

26.0 (62.3) 

r B (au) 

0.7 (1.3) 

0.4 (0.7) 

7r s i (mas) 

11.5 (244) 

62.0 (1066) 

M 

15 (51) 

31 (140) 

A (pc) 

62 (209) 

14 (44) 

D s (pc) 

360 (800) 

137 (326) 

v ( kms -1 ) 

135 (137) 

132 (134) 

cr a (pas) 

380 (570) 

81 (171) 

G magnitude 

rs 18 (18) 

rs 15 (16) 


Table 3. This records the median and the mean (in brackets) 
values of the parameters of the simulated events for the whole 
sample and for the subset of high quality events (those with mass 
uncertainty less than 50%). 


from the LF is fraught with difficulties (e.g., D’Antona & 
Mazzitelli 1994). Hence, for masses below 0.5M©, we assume 
that the MF is flat in our standard model. 


4.5 Velocity Distributions 

We draw components of source and lens velocities from tri- 
axial Gaussian distributions: 

f{V R ,V„Vt) = 
x 

For the bulge sources, the random velocity is or = cr<j, = 
a z = 100kms~ about a mean of 50kms“ (McGinn et 
al. 1989). For the disk sources, the random component has 
(Jr = 34kms _1 a $ = 21kms -1 o z = 18kms _1 about a 
mean velocity u m of 214kms~ 1 (Edvardsson et al. 1993). 
The assumption of triaxial Gaussians with constant semi¬ 
axes seems reasonable and is often made in studies of Galac¬ 
tic microlensing (e.g., Kiraga & Paczyriski 1994, Kerins et 
al. 2001). Nonetheless, a more exact calculation should de¬ 
rive the spatial variation in the means and dispersions of 
the velocities of the stellar populations by solving the Jeans 
equations. 


( 27 r) 3 / 2 CTRO><T z 

Vb (V» - Thn ) 2 

2 al 2o- 2 


exp(-/ fl 


V 

v z > 

2<t? ' 


(35) 


5 MONTE CARLO SIMULATIONS 

We first pick the source position, velocity and magnitude, 
which provides the astrometric accuracy. This is done by 
drawing randomly from the probability distributions implied 
by eqs. j27| ) and (|3S|), using the rejection method (e.g., Press 
et al. 1992). The astrometric optical depth © can be recast 
to provide the joint probability of drawing the lens distance 
and transverse velocity for a given source. We draw the mass 
of the lens from the mass function. Finally, we pick the im¬ 
pact parameter u o and the phase to by drawing uniformly 
from the rectangular area in the lens plane of size Tnt e v and 
2u a RE (illustrated in Figure 2 of Dominik & Sahu (2000)). 
This provides us with a single simulated microlensing event, 
which has a variation of the centroid shift at least as large 


as 5\/2cr a and which has a maximum during the lifetime of 
the GAIA mission. 


5.1 Distributions of the Events 

We create a set of 25000 such microlensing events, as a syn¬ 
thetic GAIA database. For each member of the sample, we 
compute the error in the estimated mass of the lens using 
eq (^5|) , which of course depends on the uncertainties in the 
relative parallax and the angular Einstein radius. From this, 
we build up a subset of 2531 events for which the error in the 
inferred mass is less than 50%. We refer to this as the sample 
of high quality events. Figures [ 12 ] and [p] show the distribu¬ 
tions of the microlensing parameters for the entire sample 
and for the subset of high quality events. The medians and 
means of the distributions are given in Table [^. Note that 
the distributions have long tails caused by high impact pa¬ 
rameter events and so we typically use medians rather than 
means, in accord with normal statistical practice. 

There are a number of interesting things to notice. The 
mean source distance is ~ 800 pc for the whole sample, and 
about 140 pc for the high quality events. The mean lens 
distances are still smaller. So, the astrometric microlensing 
signal seen by GAIA is overwhelmingly dominated by the 
very local stars. Of course, more distant microlensing events 
take place, but they can be much less commonly identified 
and still less commonly analysed to extract the parameters 
for the event. 

For the whole sample, the optimum direction is towards 
the Galactic Centre. This is not surprising, as the source 
and lens density is highest in this direction. However, the 
high quality events are found preferentially at longitudes 
l ~ 70° and 260°. These events have longer durations than 
usual. The largest component of the transverse velocity is 
the circular motion about the centre of the Galaxy. So, it 
helps to be looking in directions in which the circular ve¬ 
locity is almost entirely along the line of sight. Our model 
of the Galaxy is axisymmetric, so the peaks in the dashed 
histogram of longitude are roughly symmetric about i = 0°. 
A similar effect can be seen in the distribution of Galactic 
latitude. Here, the high quality events have a broader dis¬ 
tribution than the whole sample. A longer duration event 
is favoured by a somewhat higher latitude as the velocity 
component transverse to the line of sight is diminished. 

In understanding the distribution of accuracies cr a , we 
recall that the events in the whole sample have been chosen 
to have at least a 5\/2<r a variation in the centroid shift. If 
the source is faint, then the accuracy of the astrometry is 
poor and very few events can satisfy this stringent criterion. 
We recall from eq (113) that <r a is proportional to #e and 
varies inversely with u 2 . So, the high quality sample is still 
more strongly biased to nearby luminous stars for which er a 
is small. 

There are two ways to obtain a high quality event. Ei¬ 
ther the astrometric accuracy is good or the microlensing 
event has optimum characteristics, namely that the angular 
Einstein radius is large, the duration is long and the pro¬ 
jected Einstein radius ?e is comparable to the size of the 
Earth’s orbit around the Sun. All these characteristics are 
evident on comparing the full and the dashed histograms in 
Figure [b| The median duration of the high quality events is 
~ 350 days, nearly a factor 10 times larger than the whole 
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Source distance, D, 




Source galactic longitude, 


degre 

Transversal velocity 




Source galactic latitude, b 



degre 

Astrometric accuracy 



Figure 12. Distributions of source distance D s , lens distance D\. source coordinates (t. b), transverse velocity v and astrometric accuracy 
o a for the simulated datasets. The histogram in full lines refers to the whole sample, dotted lines to the sub-sample of high quality events. 
The distributions are all normalised to a maximum value of unity. The medians are indicated by vertical bars. 


sample. The projected Einstein radius ?e already has a me¬ 
dian of 0.7 au for the whole sample, and it is still smaller at 
0.4 au for the high quality events. The median angular Ein¬ 
stein radius is nearly a factor of 4 larger for the high quality 
events. Note that the median impact parameter for the high 
quality events is larger than that for the whole sample. This 
is because uo ~ G #E/oa. For a high quality event, it helps 
to have large angular Einstein radius and good accuracy, so 
the distribution of impact parameters has a long tail. 

Table ^ shows the percentage of events for which 
M,tE,0E and 7r s i can be recovered to within 50%. From the 
table, we see that over 40% of all the events will have good 
estimates for Ie and 6 *e from the GAIA astrometry alone. 
Current ground-based programs like MACHO and EROS 
attempt to recover the characteristic masses of lenses from 
estimates of Ie alone for samples of a few tens to hundreds 
of events. By contrast, GAIA will provide a much larger 
dataset of ~ 10000 microlensing events with good estimates 


of both Ie and #e- The typical locations and velocities of 
the sources and lenses can be inferred for the ensemble us¬ 
ing statistical techniques based on Galactic models (in much 
the same manner as MACHO and EROS do at present). 

The number of events with an amplification A corre¬ 
sponding to 5cr p and with at least one datapoint sampled on 
the photometric lightcurve by GAIA is 1260. So, ~ 5% of 
all the 25000 events have a photometric signal recorded by 
GAIA. One datapoint is of some help as it provides corrob¬ 
orating evidence that a microlensing event has taken place. 
If we ask for at least nine datapoints (which may very well 
be clustered, so there may be only three independent mea¬ 
surements), then the number of events is 427. So, for only 
~ 2% of all the astrometric events will GAIA itself obtain 
photometry which will substantially improve the parameter 
estimation. 
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Figure 13. Distributions of impact parameter « a , event duration t ae , angular Einstein radius . projected Einstein radius fg, relative 
parallax 7r s i and lens proper motion angle <fr for the simulated datasets. The histogram in full lines refers to the whole sample, dotted lines 
to the sub-sample of high quality events. The distributions are all normalised to a maximum value of unity. The medians are indicated 
by vertical bars. 


Percentage 

< 50%) 

10 % 

N(a tE < 50%) 

37 % 

7V( ff£)E < 50%) 

46 % 

A(cr„ sl < 50%) 

13 % 


Table 4. This gives the percentage of events in the whole sample 
with an error of less than 50% in the estimate of mass M, or the 
Einstein crossing time tg, or the angular Einstein radius 11 j-: or 
the relative parallax 7r s i. 


5.2 Determination of the Local Mass Function 

GAIA is the first instrument to have the capabilities of de¬ 
tecting nearby populations of very dark objects. For exam¬ 


ple, local populations of low mass (~ 0.5 Mq) black holes 
could easily have eluded identification thus far. Cool halo 
white dwarfs are also extremely difficult to detect. Existing 
programs look for faint objects with high proper motions 
characteristic of halo populations (e.g., Oppenheimer et ah 
2001). Such programs miss stars with low space motions and 
have difficulties detecting stars with extremely high proper 
motions (depending on the epoch difference of the plates). 
So, they can only set lower limits to the local halo white 
dwarf density. Samples of disk white dwarfs are again se¬ 
lected kinematically and are subject to similar biases (Knox, 
Hawkins & Hambly 1999). There may remain undetected a 
larger population of even fainter and cooler white dwarfs, as 
it is unclear whether existing surveys are probing the very 
faintest luminosities possible. Similarly, there may be an ex¬ 
tensive local population of very cool brown dwarfs. The 2 
Micron All Sky Survey (2MASS) and the Sloan Digital Sky 
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Lens mass 



Figure 14. The distribution of masses of the lenses for the whole 
sample (full histogram), the subsample of high quality events 
(dotted histogram) and the subsample of events with some pho¬ 
tometry (dashed histogram). The bars indicate the median of the 
distributions. Notice the bias at low masses in the high quality 
sample. 

Survey (SDSS) have increased the known population of lo¬ 
cal dwarf stars with spectral types later than M to over a 
hundred (e.g., Strauss et al. 1999, Burgasser et al. 2000, 
Gizis, Kirkpatrick & Wilson 2001). Nonetheless, the coolest 
brown dwarfs will have easily eluded the grasp of these and 
all current surveys. Future projects like the Space Infrared 
Telescope Facility (SIRTF) may be able to detect very old 
brown dwarfs in the solar neighbourhood by using long inte¬ 
gration times, but a large-scale survey will be prohibitively 
costly in terms of time. However, the astrometric microlens¬ 
ing signal seen by GAIA will be sensitive to local populations 
of even the dimmest of stars and darkest of these objects. 

Here, we investigate whether we can recover the local 
mass function (MF) from the astrometric microlensing signal 
seen by GAIA. We begin by examining the selection effects 
in the subsample of high quality events, the only ones for 
which individual lens masses can be recovered to reasonable 
accuracy. Figure [i~i| shows the distribution of masses for the 
entire sample and for the high quality events. The latter are 
biased towards masses greater than 0.3 Mq. This is not too 
surprising, as it is more difficult to recover the relative par¬ 
allax if the Einstein radius is small. If instead of the high 
quality events, we consider the subsample of events with at 
least one corroborating photometry datapoint, then the bias 
is reduced. Figure |I1] already suggests that GAIA’s astro- 
metric microlensing signal will be an unbiased and powerful 
probe of the local population of stellar remnants, such as 
white dwarfs and neutron stars. However, the signal seen on 
smaller mass-scales will required modelling and correction 
for selection effects before it is unbiased. 

Reid & Hawley (2000) compute the mass function (MF) 
for nearby stars from the volume-limited 8 Parsec Sam¬ 
ple using empirical mass-luminosity relationships. They ob¬ 
tain a power-law MF (/(AL) oc Af _1 ) over the range 0.1 to 
I.OMq. This is in rough agreement with the MF as inferred 
from studies of the LF of red disk stars seen with Hubble 
Space Telescope (Gould, Bahcall, Flynn 1997), which is a 
broken power-law. Above ~ 0.6Afg, the MF is nearly of the 
classical Salpeter form with f(M) oc M~ 2 ' 21 . Below 0.6 M@, 
it is flatter than Salpeter with f(M) oc AL -0 ' 9 (after cor¬ 
rection for binaries). Of course, the behaviour of the MF 


Moss Function : Flat 




Mass Function : Fall 



Figure 15. The recovery of the flat, rising and falling mass func¬ 
tions from the subsamples of high quality astrometric microlens¬ 
ing events generated from simulations. The vertical error bar as¬ 
sociated with each bin is proportional to the square root of events 
in the bin. The horizontal error bar is the size of the bin. 


becomes more uncertain as the hydrogen-burning limit is 
approached. 

We use the three MFs shown in Figure [TT] as spanning 
a range of reasonable possibilities. Above 0.5Mq, the MF 
is always derived from the Reid-Hawley luminosity func¬ 
tion. Below 0.5Mq, there are three possibilities. It may be 
flat ( f(M ) oc A/ -1 ), rising (/(A/) oc M -1 ' 44 ) or falling 
(/(A/) oc AL 0 ' 05 ). These choices are motivated by the Reid- 
Hawley LF as shown in Figure [!(]. For the falling MF, we 
choose the power-law index to be that implied by the decline 
in the LF below Mg ~ 13. For the rising MF, we choose the 
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power-law index by assuming that the rise at Mg « 11 con¬ 
tinues to fainter magnitudes and that these stars are miss¬ 
ing from the data because of incompleteness. Of course, the 
flat MF is so called because it has equal numbers of ob¬ 
jects in each decade of mass and so appears flat in a plot 
of log TV lf(M) versus logM. For each of the three MFs, we 
generate samples of 25000 astrometric microlensing events 
using Monte Carlo simulations. We extract the high quality 
events and compute the mass uncertainty using the covari¬ 
ance analysis. (In actual practice, the high quality events 
would be selected on the basis of the goodness of their \ 2 
fits). We build up the MF as a histogram. The vertical error 
bar associated with each bin is proportional to the square 
root of events in the bin. The horizontal error bar is the size 
of the bin. 

The three cases are shown in Figure |l^ with the solid 
line representing the underlying MF. The simulated data- 
points with error bars show the MFs reconstructed from the 
high quality events. It is evident that GAIA can easily dis¬ 
tinguish between the flat, rising and falling MFs. The MFs 
are reproduced accurately above ~ 0.3 Mq . Below this value, 
the reconstructed MFs fall below the true curves, as a conse¬ 
quence of the bias against smaller Einstein radii. However, 
this does not compromise GAIA’s ability to discriminate 
between the three possibilities. In practice, of course, simu¬ 
lations could be used to re-calibrate the derived MFs at low 
masses and correct for the bias. 

We have also carried out simulations with MFs con¬ 
taining spikes of compact objects, such as populations of 
~ 0.5Mq white dwarfs. They lie in the mass regime to which 
GAIA’s astrometric microlensing signal is most sensitive. So, 
such spikes stand out very clearly in the reconstructed MFs. 
We conclude that one of the major scientific contributions 
that GAIA can make is to determine the local MF. Mi¬ 
crolensing provides the only way of measuring the masses 
of individual objects irrespective of their luminosity. GAIA 
is quite simply the best survey instrument to carry out an 
inventory of masses in the solar neighbourhood. 


6 STRATEGY 

6.1 Detection of Events 

All sources detected by GAIA will be fit using the standard 
algorithms for the 6 astrometric parameters, namely 0 s ,o, fJ-s 
and P s . This is described in, for example, volumes 1 and 3 
of The Hipparcos and Tycho Catalogues (ESA 1997). So, the 
first cut on the dataset is to ask for large residuals in the 
standard astrometric parameters. Of course, microlensing is 
just one cause of large residuals out of many. Other possi¬ 
bilities include multiple companions, eclipsing systems, the 
mottled or spotted surfaces of stars and variability-induced 
movers or VIMs (Wielen 1996). VIMs were first discovered 
in Hipparcos data, as they give spurious negative parallaxes 
when fit to the standard astrometric model. In VIMs, it is 
hypothesised that the photometric centre of a close double 
star moves on the sky as one of the components has variable 
luminosity. 

Suppose we detect the event photometrically as well as 
astrometrically. Then, GAIA’s spectroscopy and multicolour 
photometry can help us distinguish microlensing events from 




Figure 16. This shows contours of “microlensing noise” in as¬ 
trometry (upper panel) and photometry (lower panel). The num¬ 
ber within each contour refers to the total number of stars inside 
the shaded region enclosed between the contours that suffer an 
astrometric distortion or photometric brightening exceeding the 
mission target accuracy. The numbers are in units of 10 s in the 
upper panel and 10 2 in the lower panel. 


Astrometric, photometric distortion 



Figure 17. This shows the distribution of deviations due to mi¬ 
crolensing at the midpoint of the mission in units of the astromet¬ 
ric and photometric mission target accuracy. The distributions are 
all normalised to a maximum value of unity. (For agiven source, 
the mission target accuracy is that achieved at the end of the 
GAIA mission. For astrometry, this is better than the instanta¬ 
neous accuracy listed in Table M by a factor of ~ 10). 
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other possible causes of astrometric shifts. GAIA’s spec¬ 
troscopy is reasonably accurate for objects brighter than 
V ~ 15, although it covers a short wavelength range near the 
Ca triplet. For objects fainter than V ~ 15, GAIA’s muli- 
colour photometry will be more useful. If the lens is dark, 
then the colours before, during and after the microlensing 
event are the same. Such is not true of spotted stars, VIMs 
and resolved binaries, which will show colour changes over 
the course of time. If the lens is luminous, then the colour 
does change over the course of the microlensing event. How¬ 
ever, GAIA will detect such a luminous lens and so the lens 
proper motion and parallax will be available from a stan¬ 
dard astrometric fit. This means that we can correct for 
the gravitational deflection of the source directly from the 
data. The multicolour photometry can also help us reduce 
the false detection rate. An unblended microlensing event is 
achromatic and does not repeat, whereas most other types 
of variability are chromatic and do repeat. Even long-period 
variables have periods of at most ~ 600 days (Cox 2000) and 
thus are short enough for GAIA to recognise as repeating. 

In propitious circumstances, close companions can be 
distinguished from microlensing by analysis of the rate of 
change of the astrometric excursion with time. At times far 
from the closest approach, the source motion is slow in a mi¬ 
crolensing event, whereas at or near the closest approach, the 
source rapidly traverses the perturbation ellipse. Although 
a nearby companion also causes an elliptical distortion to 
be superposed on the proper and parallactic motion, the 
traversal is harmonic in time and repeats. A close compan¬ 
ion will also induce radial velocity variations, which may 
be detectable with GAIA’s spectroscopy. If the stars in the 
binary are of similar luminosity but different colour, then 
there may be color variations during the centroid motion as 
well. 

If the closest approach, which corresponds to the max¬ 
imum astrometric deviation, does not occur during the life¬ 
time of the GAIA mission, then it will be difficult to be cer¬ 
tain that the centroid shift is caused by microlensing. This 
is why, in this paper, we have restricted attention only to 
events which satisfy the stringent criteria of a 5y/2a a varia¬ 
tion in the centroid shift and a maximum during the mission 
lifetime. 


6.2 Noise caused by Microlensing 

“Microlensing noise” is the name given to the weak effects of 
typically large impact parameter microlensing events. This 
can cause fluctuations in the positions of stars. There is some 
analogy with the effects of turbulence in the upper atmo¬ 
sphere, which causes wavefronts to become corrugated and 
hence causes sources to blur and twinkle. Microlensing noise 
is a fundamental limit to the accuracy to which we can know 
the positions of sources. For example, Sazhin et al. (1998, 
2001) have pointed out that fluctuations of the angular posi¬ 
tions of reference extragalactic radio and optical sources will 
be caused by the microlensing effects of stars in the Galaxy. 
Such angular fluctuations may range up to hundreds of mi- 
croarcseconds. Sazhin et al. (1998) found a good example of 
a quasar with large positional uncertainties caused by the 
lensing effects of a nearby Hipparcos star; the distance of the 
star is very close, just 50 pcs. 

Figure [b] shows a contour plot of number of sources 


which have an astrometric or photometric variation caused 
by microlensing of the same order of magnitude as the mis¬ 
sion target accuracy. (Let us recall the mission target accu¬ 
racy is that achieved at the end of the mission and so is the 
instantaneous astrometric accuracy given in the final line of 
Table [I] divided by ~ 10). The total number of sources hav¬ 
ing an astrometric deviation exceeding the mission target 
accuracy is 2.2 x 10 5 . The positional measurement of one 
source in every twenty thousand is affected by microlensing 
noise at any instant. This number is obtained by computing 
(2.2 x 10 s /10 9 ) x (f ae }/5 years. Many of these events will be 
long duration events. Any astrometric deviation will often 
be only slowly varying over the mission lifetime and so this 
will not compromise the proper motion and parallax tar¬ 
get accuracies for GAIA. For photometry, the instantaneous 
number of stars affected by microlensing noise is ~ 1.8 x 10 3 . 
Figure 0 shows the distribution of these disturbances in 
units of the mission target accuracy. To be included on this 
Figure, an event must already have a deviation greater than 
the mission target accuracy. It is reassuring that the mean of 
the distribution of photometric deviations is strongly peaked 
near unity, so the blurring in photometry is mostly insignif¬ 
icant. 

If the distribution of errors caused by instrumental ef¬ 
fects is Gaussian, then the number of events with error less 
than 0. In a is just ~ 10%. With the more realistic assump¬ 
tion of a non-Gaussian error distribution with higher wings 
the number of such deviation will be still smaller. So, if we 
assume that microlensing noise is a random error source, 
then its effects for GAIA are negligible, as in ~ 90% of cases 
the microlensing noise is less than that expected from in¬ 
strumental sources. 

Another way in which the microlensing noise manifests 
itself is in distortion of the point spread function (PSF). 
GAIA’s PSF may be approximated as a biaxial Gaussian 
with FWHM of 300 mas and 150 mas. To investigate the 
likely size of the errors, we convolve the PSF with the two 
micro-images corresponding to an event with Einstein ra¬ 
dius 9e = 5 mas. These micro-images are taken as pure 
5-functions. The convolved PSF is then fit with the stan¬ 
dard biaxial Gaussian. This procedure of course carries an 
error, which we estimate by computing the area difference 
between the true PSF and fit. This is divided by the area 
under the PSF to give the fractional error caused by mi¬ 
crolensing. Table list the percentage error as a function 
of impact parameter. Although the error is small, it will be 
important as GAIA is aiming for a photometric accuracy of 
fractions of a percent. 


6.3 Photometry and Astrometry 

Figure [l8| shows the uncertainty in the relative parallax as a 
function of astrometric accuracy for events with astrometry 
alone and events with astrometry and photometry. The lat¬ 
ter curves are calculated using the covariance analysis, but 
now incorporating the photometric measurements as well 


-2 d6 a (tk) dd a {tk) -2 
-fa — d^- + aA 

k =1 


dA(tk) dA(tk) 

dai da,j 


where A is the photometric amplification and 
a A = 10°- 4<Tp - 1. 


(36) 


( 37 ) 
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u 0 0.25 0.5 1 5 

% error 0.28 0.26 0.22 0.14 


Table 5. This table lists the percentage error in the estimated point spread function against the impact parameter of the microlensing 
event. 


Parallax uncertainty vs. accuracy 



Figure 18. This shows the percentage error in the parallax of 
the lens as a function of the astrometric accuracy without and 
with additional GAIA photometry. The parameters for this event 
are as in Figure 1, but the impact parameter u is 2.5. 

The figure makes the point that just by using GAIA’s pho¬ 
tometry, it is possible to improve the accuracy of the param¬ 
eter estimation by a factor of ~ 10. This however will only 
be viable for the ~ 3% of events for which GAIA records 
a number of datapoints on the photometric lightcurve (see 
Section 5.1). 

Can GAIA trigger detection of microlensing events 
which can subsequently be followed-up on with ground- 
based telescopes? This is a real challenge. GAIA provides 
a one-dimensional stream of astrometric data, and so must 
operate for roughly two years to gain enough accuracy on 
parallactic motions. After two years, preliminary astromet¬ 
ric fits may be constructed, but their accuracy will be lim¬ 
ited. The best opportunities for detecting events will come 
in the latter part of the mission, when the iterative astro- 
metric solution becomes more accurate. Hence, there is a 
need for at least one dedicated lm class telescope (or bet¬ 
ter still a world-wide network of such telescopes) to monitor 
fields in the bulge and spiral arms so that the datasets can 
be correlated with GAIA’s database. The error is reduced 
by typically better than a factor of 10 when astrometric 
data is supplemented by densely-sampled lightcurves from 
the ground (see e.g., Gould & Salim 1999). 

6.4 The FAME mission 

A related satellite is the Full-Sky Astrometric Mapping Ex¬ 
plorer or FAME [|], which is scheduled for launch in 2004 and 
has a mission lifetime of 2.5 years. It is based on the same 
principles as the Hipparcos satellite. It is less ambitious than 
GAIA, as it monitors mainly sources within about 2.5 kpc of 
the Sun. FAME will provide astrometry on all objects with 

t http://www.usno.navy.mil/FAME/ 


6 < V < 16. For a source with V = 10, the accuracy of 
individual measurements cr a is 750 p,as and the mission tar¬ 
get accuracy is 36 /ras. By comparing with Table [jj we see 
that, roughly speaking, the accuracy of FAME’S individual 
measurements is a factor of 10 worse than that of GAIA. 

The all-sky averaged astrometric optical seen by FAME 
is ~ 2 x 10 -6 . There are ~ 10 8 stars with 6 < V < 16, 
which means that FAME will measure ~ 200 astrometric 
microlensing events (with the characteristics that they have 
a 5y/2 <r a deviation and that the maximum occurs during the 
mission lifetime). This number is low, because the mission 
lifetime is short and the individual measurements are noisy. 
It is already clear that the high quality sample will be too 
small to deduce any useful information about the local mass 
function. 


7 CONCLUSIONS 

GAIA is an ambitious astrometric survey satellite that is 
planned for launch no later than 2012 by the European Space 
Agency (see ESA 2000; Perryman et al. 2001). GAIA scans 
the whole sky, performing multi-colour and multi-epoch pho¬ 
tometry and spectroscopy, as well as astrometry. The refer¬ 
ence frame is tied to a global astrometric reference frame 
defined by extragalactic objects. For bright sources, GAIA 
has the capabilities to perform individual astrometric mea¬ 
surements with microarcsecond accuracy. Even for objects 
as faint as V « 20, the positions, proper motions and par¬ 
allaxes will be recovered to within at worst a few hundred 
microarcsecs at the end of the mission. 

GAIA can measure microlensing by the small excursions 
of the light centroid that occur during microlensing. We use 
a stringent definition of an astrometric microlensing event 
as one with a significant variation (5\/2cr a ) of the centroid 
shift, together with a closest approach during the lifetime 
of the GAIA mission. The all-sky averaged astrometric mi¬ 
crolensing optical depth is ~ 2.5 x 10 -5 . This means that 
~ 25000 sources will exhibit astrometric microlensing events 
with such characteristics during the course of the mission. 
GAIA can also measure microlensing by the photometric 
brightening that accompanies low impact parameter events. 
The all-sky averaged photometric optical depth is ~ 10 - ', 
so there are ~ 3600 photometric microlensing events during 
the five year mission lifetime, most of which are undetectable 
because of the poor sampling. Consequently, very few astro- 
metric events are measured photometrically as well. In fact, 
only ~ 1260 astrometric events have at least one datapoint 
sampled on the photometric lightcurve. Only for ~ 427 of 
the astrometric events will GAIA itself obtain enough pho¬ 
tometric datapoints to improve substantially the characteri¬ 
sation of the event. Let us again emphasise that these num¬ 
bers only refer to numbers of microlensing events with spec¬ 
ified characteristics. Any statement as to the actual number 
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of microlensing events that GAIA will find must take into 
account the number of false detections thrown up by the 
identification algorithm. 

Our event statistics refer only to the signal provided 
by stellar lenses in the Galactic disk. We have not taken 
into account either stellar lenses in the Galactic bulge or 
dark objects in the halo (popularly called Machos). This is 
reasonable, as the cross-section for astrometric microlensing 
favours nearby lenses. The most valuable events are those for 
which the Einstein crossing time £e, the angular Einstein ra¬ 
dius 9e and the relative parallax of the source with respect 
to the lens 7T s i can all be inferred from GAIA’s datastream. 
The mass of the lens then follows directly. If the source dis¬ 
tance is known - for example, if GAIA itself measures the 
source parallax - then a complete solution of the microlens¬ 
ing parameters is available. Of these quantities, it is the 
relative parallax that is the hardest to obtain accurately. A 
covariance analysis is used to follow the propagation of er¬ 
rors and establish the conditions for recovery of the relative 
parallax. This happens if the angular Einstein radius 9e is 
large and the Einstein radius projected onto the observer 
plane ve ~ 1 au so that the Earth’s motion about the Sun 
gives a substantial distortion. It is also aided if the source is 
bright so that GAIA’s astrometric accuracy is high and if the 
duration of the astrometric event is long so that GAIA has 
time to sample it fully. These conditions favour still further 
lensing populations that are very close. 

Monte Carlo simulations are used to establish the char¬ 
acteristics of the ~ 10% of events for which GAIA can re¬ 
cover the mass of the lens to good accuracy. The typical 
lens distance is ~ 50 pc and the typical source distance is 
~ 300 pc. The highest quality events seen by GAIA are 
overwhelmingly dominated by very local lenses. 

We conclude that one of the major scientific contribu¬ 
tions of microlensing studies with GAIA will be the deter¬ 
mination of the mass function in the solar neighbourhood. 
Of course, direct mass measurements are presently possible 
just for binary stars with well-determined orbits. Microlens¬ 
ing is the only technique which can measure the masses of 
individual stars. GAIA is the first instrument with the abil¬ 
ity to survey the astrometric microlensing signal provided 
by nearby lenses. We have used Monte Carlo simulations to 
show that GAIA can reconstruct the mass function in the 
solar neighbourhood from the sample of its highest quality 
events. This works particularly well for masses exceeding 
~ O.SMq. Below 0.3Mq, the reconstructed mass function 
tends to underestimate the numbers of objects, as the high¬ 
est quality events are biased towards larger angular Einstein 
radii. 

If there are local populations of low mass black holes, 
or very cool halo and disk white dwarfs or very old brown 
dwarfs, then they will have easily eluded detection with 
available technology. However, the astrometric microlensing 
signal seen by GAIA will be sensitive to local populations 
of even the dimmest of these stars and the darkest of these 
objects. GAIA is the first instrument that has the potential 
to map out and survey our darkest neighbours. 
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